Glass transition in models with controlled frustration 
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A class of models with self-generated disorder and controlled frustration is studied. Between 
the trivial case, where frustration is not present at all, and the limit case, where frustration is 
present over every length scale, a region with local frustration is found where glassy dynamics 
appears. We suggest that in this region, the mean field model might undergo a p-spin like transition, 
and increasing the range of frustration, a crossover from a 1-step replica symmetry breaking to a 
continuous one might be observed. 
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The study of spin glasses and glass formers has shown 
that these systems present similar complex dynamical 
behaviors. In particular the dynamical equations of Jo- 
spin models 0], a generalization of the spin glass model 
in mean field, coincide above the dynamical transition 
temperature with the mode coupling equations for su- 
percooled liquids ^2J. However the connection in finite 
dimension between spin glasses and glass formers is not 
completely clear. On one hand, spin glasses undergo a 
thermodynamic transition at a well defined temperature, 
where the nonlinear susceptibility diverges. The systems 
that show a transition of this kind, in spite of very dif- 
ferent microscopic structures, have two essential common 
characteristics: The presence of competitive interactions 
(frustration) and quenched disorder. On the other hand, 
glass formers are a class of systems where disorder is 
not originated by some fixed external variables, but is 
"self-generated" by the particle positions and orienta- 
tions. Moreover there is no sharp thermodynamic transi- 
tion characterized by the divergence of a thermodynamic 
quantity analogous to the nonlinear susceptibility. In or- 
der to clarify the connection and the differences between 
glasses and spin glasses, and to investigate the roles of 
disorder and frustration in the behaviors observed, in the 
present paper we study a class of models with annealed 
interactions and controlled frustration. 

We consider a diluted spin glass, the frustrated lat- 
tice gas (FLG) _3|, constituted by diffusing particles, and 
therefore suited to study quantities like the diffusion co- 
efficient, or the density autocorrelation functions, impor- 
tant in the study of liquids. The Hamiltonian of the 
model is: 



(1) 



where Si 
variables 



= ±1 are Ising spins, = 0, 1 are occupation 
and eij = ±1 are ferromagnetic and antifer- 
romagnetic interactions between nearest neighbor spins. 
This model was studied both for quenched 0, 01 and 
annealed interactions "s"] : In the quenched case the inter- 
actions, eij, are quenched variables randomly distributed 
with equal probability; in the annealed case evolve in 
time. 



In the limit fj,/T goes to infinity all sites are occupied 
(ui = 1 for each site i), and the quenched model repro- 
duces the Ising spin glass. In the other limit, T/J = 0, 
the model, eq. has properties recalling a "frustrated" 
liquid. Indeed the first term of the Hamiltonian implies 
that two nearest neighbor sites can be freely occupied 
only if their spin variables satisfy the interaction, that 
is if eijSiSj ~ 1, otherwise they feel a strong repulsion. 
Since in a frustrated loop the spins cannot satisfy all 
the interactions, in this model particle configurations in 
which a frustrated loop is fully occupied are not allowed. 
The frustrated loops in the model are the same of the 
spin glass model and correspond in the liquid to those 
loops which, due to geometrical hindrance, cannot be 
fully occupied by the particles. 

Here we study a class of annealed FLG models defined 
by the following partition function: 
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where H. is given by eq. (^3), the sum X^jo-} over all 
the possible configurations of spin and particles {it} = 
{Si, Tii}, and the sum X]{e} over all the possible inter- 
action configurations such that the annealed averages of 
the frustrated loop number coincide with the quenched 
ones for every length of the loops until a maximum value, 
0, 4, 6, 8, ... on a cubic lattice). Varying 
fmax from zero to infinity a class of models with con- 
trolled frustration is obtained. In the limit, rmoK — 0, 
frustration is not present at all: the model Q is equiv- 
alent to a lattice gas with a repulsion between nearest 
neighbors, and without frustration and correlations be- 
tween spin, and no thermodynamic transition is present. 
In the other limit, rmax goes to infinity, frustration is 
present over every length scale, as in the quenched case: 
We expect that the partition function, eq. (0), coincides 
with the quenched one, namely the model undergoes a 
spin glass- like transition j^, as shown in Ref. €] for the 
Ising spin glass model 1^ . For intermediate values of rmax 
a class of models with local frustration and self-generated 
disorder is obtained. 

In the present paper the models, eq. have been 
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studied for rmax = 4 and Vmax = 6: In the first case 
a trivial dynamical behavior is observed, with one step 
relaxation functions and a smooth increasing of relax- 
ation time as function of density; in the second case a 
dynamical behavior very similar to that of glass formers 
is instead found. We therefore observe that by increas- 
ing the degree of frustration the system moves from a 
liquid-like to a glassy-like behavior. 

The system was simulated using Monte Carlo tech- 
niques over a cubic lattice of linear size L = 8. At the 
beginning the interactions, e.y = ±1, arc randomly dis- 
tributed with equal probability. At each step of the dy- 
namics an attempt to move a particle to a nearest neigh- 
bor site (the spin is flipped with a probability equal to 
1/2) is alternated with an attempt to exchange two near- 
est neighbor interactions: in the limit here considered, 
T/ J = 0, a particle can move to a nearest neighbor site 
only if its spin satisfies the interactions with all the new 
nearest neighbors, and an interaction can be changed 
only if at least one of its extreme is empty. The frus- 
trated loop numbers of any fixed length until r^ax are 
independently kept constant during the dynamics. 

At a given value of the density, p, we calculated 
the two-time relaxation function of the self-overlap, 

i{tw)Si{t)ni{t), where Np 
is the number of particles, and the average is done 
over 8 — 32 different realizations of the system. For values 
of long enough, the system reaches a stationary state, 
where the time translation invariance is recovered, i.e., 
C{t,tyj) = C{t — tm)- In this time region we calculated 
the equilibrium relaxation function of the self-overlap: 

W = ■^J2^S,{t)n,{t)S,it + t')n,{t + t')), (3) 

and the dynamical nonlinear susceptibility 0: 

Xit)=Npi{q^)it)-{qy{t)). (4) 

Here (• ■ •) is the time average on time t' . For each density 
the quantities of interest are averaged over 8 — 32 different 
realizations of the system, and the errors are calculated 
as the fluctuations over this statistical ensemble. 

We first consider the model with r^ax = 4: Since r^ax 
equals the loop minimum length, the interactions evolve 
under the constraint that the number of frustrated loops 
of length 4 is kept constant. In this case, both {q){t) and 
x{t) show a liquid-like behavior also at high densities: 
{q){t), plotted in Fig. relaxes with a one step decay 
well fitted by a stretched exponential function, f{t) — 

A exp|— (r)*^!: with S ~ 0.93 not depending on the 

density; and xi^) tends to a plateau, which smoothly 
increases as a function of density. We suggest that this 
behavior might be due to the fact that there is only a few 
residual frustration on loops of length greater than 4, and 
frustration is too local to originate a slow dynamics. 

A very different behavior is instead observed in the 
model with Vmax = 6 (where the interactions evolve un- 
der the constraint that the number of frustrated loops of 



<q>(t) 




FIG. 1: The relaxation functions of the self-overlap, {q){t), for 
r„„^ = 4 at the densities p = 0.63, 0.70, 0.74, 0.76 (from 
left to right). The continuous curves in figure are stretched 
exponential functions with S = 0.93. 
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FIG. 2; The relaxation functions of the self- 
overlap, {q){'t), for Tmax = 6 at the densities 
p = 0.48, 0.54, 0.59, 0.63, 0.66, 0.71, 0.73 (from left 
to right). 



length 4 and 6 are independently kept constant.) In Fig. 
12 (q) (t) are plotted at different values of the density. At 
high density two step decays appear and the curves are 
well fitted at intermediate times by the forms predicted 
near the dynamical transition by the mode coupling the- 
ory (MCT) HH (see Fig. ^. In Fig. H the scaled 
relaxation functions of the self-overlap, {q){t/T)/fQ, are 
plotted as function of the scaled time, t/r. At long times 
the curves collapse onto a single master function, well fit- 
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FIG. 3: The relaxation functions of the self-overlap, {q){t), 
for Tmax = 6 at the density p = 0.69. The curve in figure 
is the /3— correlator of the MOT with exponent parameters 
a = 0.327 and b = 0.641, and plateau height /o = 0.84. 



FIG. 5: The relaxation time, r, and the inverse diffusion co- 
efficient, , (the empty circles in figure) are plotted as 
function of pr — p. The curves are the power laws obtained 
as fitting functions. 
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FIG. 4: The scaled relaxation functions of the self-overlap, 
(g)(t/r)//o, as a function of the scaled time, t/r, at the den- 
sities p = 0.54, 0.59, 0.63, 0.66, 0.71 (from right to left). 
The continuous curve is a stretched exponential function with 
S = 0.71. 



ted by a stretched exponential (the continuous curve in 
figure), and the relaxation time, r, diverges as a power 
law, (pr-p)"'^^ with pr ~ 0.79 ±0.02 and 7^ = 4.9 ±0.7 
(see Fig. 01 . 

The dynamical nonlinear susceptibility, x(i), shown in 
Fig. El presents a maximum at a time, t* , which we inter- 
pret as the relaxation time of the interactions: Until t* 
the dynamical nonlinear susceptibility increases as if the 



FIG. 6: The dynamical nonlinear susceptibility, x(^)i for 
rmax = 6 at the densities p = 0.48, 0.54, 0.59, 0.69 (from 
left to right). 



environment were quenched, and only for t > t* the inter- 
actions evolve and x{t) can decrease to the equilibrium 
value. A dynamical nonlinear susceptibility with a max- 
imum is typical of glassy systems [Sl llOj . As in Lennard- 
Jones liquids we found that the value of the maximum, 
x(t*), which diverges in the p-spin model as the dy- 
namical transition is approached from above, has instead 
a maximum: We suggest that in the present case this 
behavior might due to the fact that the system becomes 
less and less frustrated as the density increases [Tll[ . 
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Finally we calculated the particle mean square dis- 

placement, {Ar^)it) = l/TV^ ^f" ((r,(t + t') - r.(i'))'), 
where ri{t) is the position of the i— th particle at the time 
t. At high density the mean square displacement is not 
linear at short time, and the diffusion coefhcient, Z?, is 
calculated from the long time regime of the mean square 
displacement via the relation, D = limt^oo(Ar^)(i)/6t. 
The diffusion coefficient, shown in Fig. [Sj is well fitted 
by a power law, {po — p)'*" , with po — 0.80 ± 0.01 and 
7d = 2.46 ± 0.19. The critical density, po, obtained 
in this way coincides with the value, p^ = 0.79 ± 0.02, 
where the relaxation time, r, diverges; the exponent, 7c 
is instead not consistent with = 4.9 ± 0.7. 

In conclusions the properties of the annealed models, 
eq. Q, strongly depend on the value of rmax- In par- 
ticular by increasing the degree of frustration a crossover 
from a liquid-like to a glassy-like behavior is found. We 
suggest that the model with rmax = 6, where a dynamical 
glass transition is observed, in mean field might undergo 



a p-spin like transition with a 1-step replica symmetry 
breaking in the spin overlap distribution. Moreover it is 
reasonable to expect that, by further increasing the de- 
gree of frustration, the annealed partition function, eq. 
(0), might tend to the quenched one, where a spin glass- 
like transition and the development of a continuous 
replica symmetry breaking in the spin overlap distribu- 
tion is found, and a crossover from glassy-like to spin 
glass-like behavior might be observed. 
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